Astronomy & Astrophysics manuscript no. 12666shu 


© ESO 2009 


September 1, 2009 





o 
o 

(N 

Oh 

c/3 



X 

Oh 

6 

o3 



> 

in 

o 
q 

ON 

o 

o> 
o 



X 



INTEGRAL and Swift/XRT observations of IGR J19405-3016 

(Research Note) 

Shu Zhang 1 , Yu-Peng Chen 1 , Diego F. Torres 2 , Jian-Min Wang 1,3 , Ti-Pei Li 1,4 , Jun-Qiang Ge 1 

1 Key Laboratory for Particle Astrophysics, Institute of High Energy Physics, Beijing 100049, China 

2 ICREA & Institut de Ciencies de l'Espai (IEEC-CSIC), Campus UAB, Facultat de Ciencies, Torre C5-parell, 2a 
planta, 08193 Barcelona, Spain 

3 Theoretical Physics Center for Science Facilities (TPCSF), CAS 

4 Center for Astrophysics, Tsinghua University, Beijing 100084, China 

Received / Accepted 

ABSTRACT 

Aims. IGR J19405-3016 is reported in the 3rd IBIS catalog as one of its lowest significance sources (~ 4.6 a under an 
exposure of about 371 ks). This leads to a caveat in multi- wavelength studies, although the source was identified in 
the optical as a Seyfert 1. The currently available INTEGRAL data on the source have increased to an exposure time 
of ~ 1400 ks, which stimulates us to investigate the reality of this source again by using all the available data from 
INTEGRAL and Swift/XRT. 

Methods. We analyzed all available observations carried out by the International Gamma- Ray Astrophysics Laboratory 
(INTEGRAL) on IGR J19405-3016. The data were processed by using the latest version OSA 7.0. In addition, we 
analyzed all the available Swift/XRT data on this source. 

Results. We find that IGR J19405-3016 has a detection significance of ~ 9.4 a in the 20-60 keV band during the 
observational period between March 2003 and March 2008. This confirms a real source detection reported previously. 
The source position and error location are therefore updated. The source is found to be constant over years at the hard 
X-rays. We analyzed the Swift/XRT observations on IGR J19405-3016 as well, and find that the spectrum can be fitted 
with a simple power law model. Over the three XRT observations, the source flux varied by up to 39% from the average, 
and the spectrum is generally soft. The combined XRT/ISGRI spectrum is well fitted with a simple power law model 
(photon index 2.11±0.03) with a column density fixed at 8.73xl0 20 atoms/cm 2 . Such a photon index is consistent with 
the mean value 1.98 (dispersion 0.27) obtained from Swift /BAT AGN samples at 14-195 keV. The spectral slope of 
IGR J19405-3016 is softer than the average spectral slope found elsewhere. A similar discrepancy is found with other 
results regarding Seyfert 1 AGNs. A possible explanation for this simple spectral description may be that the low level 
of the column density allows for the 'true' spectrum to appear at soft X-rays as well. 

Key words. X-rays: individual: IGR J19405-3016 



^ ' 1. Introduction IGR J19405-3016 was first reported as a low-significance 

, , ,.-,.,„, source in the 3rd IBIS/ISGRI catalog (Bird et al. 2007). 

S^/t^t t y nt y; on J e S0 T e onn a Jf I 'V ?~ The source was detected at 4.6 a level under an exposure 

£ . IBIS/ISGRI catalog (Bird et al . 2007). Among them, 171 of ^ m kg ^ ^ band 20 _ 4Q keV The bwest gi 

are Galactic accreting systems, 122 extragalactic objects, nificance for a source in Bird > s catal is 45 a There f re, 

and 113 are sources of unknown nature. Most of the un- IQR J19405 _ 3016 is a the few sources that were re _ 
classified sources were identified through optical and near- ted b INT EGRAL/ISGR1 wit h the lowest significances, 

onnf on^?^?^^^ 1 et , t'Jl ^ P ' Bird et al. (2007) point out that the catalog sources de- 

2 ? on C d 0( 2 d ' I 008 ' 2 ,°° 9 i Ch f I 6t aL 200 k ; 6t tected with significance above 5 a have a probability of 

al 2008) For the newly identified sources, about 55% are lesg than 1 nt of bei rio and detections with 

AGNs (almost equally divided into Seyfert Is and 2s) 32% Wr si mficance have a hi hcr probabmty of bein g unrea l. 

are X-ray binaries and 12% are CVs (Masetti et al. 2009). Accordin i a i tho ugh a Seyfert 1.2 was found as an opti- 

The optical identifications are performed by first searching ^ counte art around IGR ji 94 05-3016 by Masetti et al. 

m the error location ol IINTEGRAL/ISGRI for a possible (2008)j & cayeat wag daimed therein ^ ^ gense ^ ^ 

counterpart at the soft X-ray observations and, accordingly, hard x _ ray detection itself be spuriouSj and henC e the 

the error location can be largely improved to a level that i i • i ,- c . j . 

„ . , , . . . t i r i n i • i -n proposed optical identihcation may not correspond to an 

allows for looking into optical data for the hnal identinca- , > nn ,„ „ , / 1 \ /r „ „ ++ - . „i onno^ 
t ^.t^ . r , r»rmn^,n a t i, ua l counterpart ( Masct t i et al . 2008 ) . 

tion. IGR J19405-3016 is one of the INTEGRAL weakly H v ' 

identified sources studied by Masetti et al. (2008) in their After the discovery of the source by INTEGRAL, 

most recent optical identification campaign. Swift/XRT detected a source at soft X-rays within the IBIS 

error location (Landi et al. 2007a). The XRT observation 

Send offprint requests to: Shu Zhang locates the source at RA(J2000) = 19h 40m 15.15s and 
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certainty at 90% confidence level. The X-ray spectrum is 
modeled with an unabsorbed power law with a photon in- 
dex ~ 2.1 (Landi et al. 2007a). Given that an additional 
exposure will lead to a large increment in source signifi- 
cance and that the currently available INTEGRAL obser- 
vations around IGR J19405-3016 have been accumulated 
to roughly 1400 ks, it is appropriate to re-investigate the 
hard X-ray properties of IGR J19405-3016. Also, after the 
report of Landi et al. (2007a) , further Swift /XRT observa- 
tions (roughly 1.6 ks more) were available in 2008, which 
allow for further investigation on the source at soft X-rays. 
In this paper we report the results of our INTEGRAL and 
Swift /XRT analyses on IGR J19405-3016. 

2. Observations and data analysis 

INTEGRAL (Winkler et al. 2003) is a 15 keV - 10 MeV 7- 
ray mission. The main instruments are the Imager on Board 
the INTEGRAL Satellite (IBIS, 15 keV - 10 MeV; Ubertini 
et al. 2003) and the SPectrometer onboard INTEGRAL 
(SPI, 20 keV - 8 MeV; Vedrenne et al. 2003). They are 
supplemented by the Joint European X-ray Monitor ( JEM- 
X, 3-35 keV) (Lund et al. 2003) and the Optical Monitor 
Camera (OMC, V, 500-600 nm) (Mas-Hesse et al. 2003). At 
the lower energies (15 keV - 1 MeV), the CdTe array ISGRI 
(Lebrun et al. 2003) of IBIS has a better continuum sensi- 
tivity than SPI. The satellite was launched in October 2002 
into an elliptical orbit with a period of 3 days. Due to the 
coded-mask design of the detectors, the satellite normally 
operates in dithering mode, which suppresses the system- 
atic effects on spatial and temporal backgrounds. 

The INTEGRAL observations were carried out in the 
so-called individual SCience Windows (SCWs), with a 
typical time duration of about 2000 seconds each. Only 
IBIS/ISGRI public data were taken into account, because 
the source is too weak to be detected by JEMX and SPI. 
The available INTEGRAL observations, when IGR J19405- 
3016 fell into the 50% coded field of view of ISGRI (offset 
angle less than 10 degrees), comprised about 566 SCWs, 
adding up to a total exposure time of ~ 1400 ks (until 
March 22, 2008). IGR J19405-3016 therefore had roughly 
1000 ks of exposure more than used in the previous report 
(Bird et al. 2007). The details of the analyzed INTEGRAL 
observations on IGR J19405-3016, including the exposure 
and the time periods, are summarized in Table 1. Most of 
these observations were carried out in the 5x5 dithering 
mode. We subdivided the data into 4 groups according to 
the observational sequence. The other observational modes 
of the selected data are staring mode and Hexagonal mode, 
the sum of which comprise only 5.6% of the whole data. 
Data from these two modes can just as well be used to 
produce a mosaic map without introducing any systematic 
errors. This has been confirmed through our consulting the 
INTEGRAL help desk. The data reduction was performed 
by using the standard Online Science Analysis (OSA) soft- 
ware version 7.0, the latest released version. The results 
were obtained by running the pipeline from the flowchart 
to the image level and the spectrum level. The flux and the 
detection significance were derived in the mosaic map at 
the source position revised in this work. 

Swift (Gehrels et al. 2004) is a 7-ray burst explorer 
launched November 20, 2004. It carries three co-aligned de- 
tectors: the Burst Alert Telescope (BAT, Barthelmy et al. 
2005), the X-Ray Telescope (XRT, Burrows et al. 2005), 



and the Ultraviolet/Optical Telescope (UVOT, Roming et 
al. 2005). We took only Swift /XRT data into account, be- 
cause BAT data were not available. The XRT uses a grazing 
incidence Wolter I telescope to focus X-rays onto a state- 
of-the-art CCD. XRT has an effective area of 110 cm 2 , an 
FOV of 23.6 arcminutes, an angular resolution (half-power 
diameter) of 15 arcseconds, and it operates in the 0.2-10 
keV energy range, providing the possibility of extending 
the investigation on the source to soft X-rays. 

There are three Swift snapshots available for IGR 
J19405-3016, each one with an exposure over 1 ks. The ob- 
servations were carried out in the photon-counting mode, 
with exposures of 7.6 ks (ID 00036657004, on July 31, 2007), 
5.9 ks (ID 00036657005, on August 7, 2007) and 1.6 ks (ID 
00036657005, on June 7, 2008), respectively. The first two 
observations were reported in Landi et al. (2007a). See de- 
tails in Table 2. We analyzed the Swift /XRT 0.3-7 keV data 
by using the latest released analysis software, provided in 
HEAsoft version 6.4. The XRT data reduction follows those 
described in Landi et al. (2007b). Here the source events 
were extracted within a circular region of radius ~ 40 pix- 
els, centered on the source position. A radius of 20 pixels 
(corresponds to 47 arcseconds) encloses about 90 percent 
of the PSF at 1.5 keV (Capalbi et al. 2005). A larger radius 
extracts more source counts, hence improves the statistics 
in spectral fittings. The spectra were fitted with XSPEC 
vl2.3.1 (Dorman & Arnaud 2001) and the model parame- 
ters estimated at 90% confidence level. 

3. Results 

3.1. INTEGRAL 

The imaging analyses show that the best source detection, 
~ 9.4 a at 20-60 keV (Fig. [TJ, was derived in the mosaic 
map of all data. The source is not detectable at lower ener- 
gies by JEMX or at higher energies by SPI. Such a detection 
is much more significant than the previous report of a 4.5 a 
signal with an exposure of 371 ks, at 20-40 keV (Bird et al. 
2007). The source flux is about 1.40±0.15 mCrab in the 20- 
60 keV band over the period 2003-2008. Under such a high- 
significance detection at 20-60 keV, the source position is 
improved to RA/Dec (J2000) = 295.0895%30.2732°, with 
a radius of 3.4 arcminutes in error circle ( 90% confidence 
level), also improved from the previous report of 5.4 ar- 
cminutes (Bird et al. 2007). 

We also investigated the source detection in individ- 
ual observational groups as listed in Table 1, in the 20- 
60 keV band. The source is detected in Revs. 0056-0244 
(MJD=52729-53291) with a significance of 4.9 a and a flux 
of 1.5±0.3 mCrab; in Revs. 0258-0371 (MJD=53332-53672) 
with a significance of 3.5 a and a flux of 1.4±0.4 mCrab and 
Revs. 0416-0498 (MJD=53805-54050) with a significance of 
7.8 a and a flux of 1.5±0.2 mCrab. One sees that the source 
is detectable in each observational group and the flux re- 
mains constant at about 1.5 mCrab level over years at hard 
X-rays. The ISGRI data after Revs. 0498 consist of only 16 
sews (Revs. 543-664), which are about 17 ks exposure in to- 
tal. The source was not detected with this exposure, and we 
have a 2-er flux upper limit of 3.3 mCrab at 20-60 keV. If, 
during Revs. 543-664, the source kept at the same flux level 
as that of the previous revolutions, it would have only been 
'clearly detected' at 1-tr level. See details for the results in 
Table 1. The source is detected in the 4 adjacent energy 
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bands of 20-25, 25-30, 30-40, and 40-60 keV, with detec- 
tion significances of 4.0, 5.4, 6.1, and 4.2 a, respectively, by 
combining all the ISGRI data. If there is a constant source 
at hard X-rays, an exposure of 1400 ks can improve the de- 
tection significance by a factor of 1.95 with respect to the 
previous report (Bird et al. 2007), which means a detection 
significance of ~ 9 a, similar to what we obtain from the 
summed data. In short, all results derived by us indicate a 
real and steady source. 

3.2. Swift/XRT 

The Swift/XRT imaging analysis shows a single source 
detected at soft X-rays within the error location of IGR 
J19405-3016. Figure 2 is the Swift/XRT image produced by 
combining the three XRT observations to show the most ac- 
curate location of the source at soft X-rays. The position is 
obtained as RA(J2000) = 19h 40m 15.00s and Dec(J2000) 
= -30d 15m 48.6s, with an error radius of 3.5 arcseconds 
at 90% confidence level, which is consistent with the report 
in Landi et al. (2007a). Over plotted are the positions of 
IGR J19405-3016, as derived in Bird et al. (2007) and in 
the present work, and their error circles. It is obvious that 
the revised position and improved error circle are more in- 
dicative of the correlation between the soft X-ray source 
detected with XRT and the hard X-ray source IGR J19504- 
3016. 

To look into the properties of the XRT source at soft 
X-rays, we carried out the spectral analysis of the three 
XRT observations. We find that data from each of the 
three observations can be fitted by using a simple power law 
model. The data in ID 00036657004 and ID 00036657005 
need an additional component of absorption, with a column 
density derived as 5-6 x 10 20 atoms/cm 2 , but this is not 
necessary for the data from observation ID 00036657006, 
probably due to a low exposure of only ~ 1.6 ks. The 
source was clearly detected in XRT observation of 2008 
June as well. The relatively low flux level and low ex- 
posure in this observation led to larger errors in model 
parameters. The reduced \ 2 were derived with values of 
< 1.0. The column density as measured directly from 
the XRT spectra of IGR J19405-3016 is rather low: a 
value around 6xl0 20 atoms/cm 2 can be even comparable 
to the Galactic column density obtained with the web ver- 
sion (9.15x 10 20 atoms/cm 2 J] and from Dickey & Lockman 
(1990) (8.73xl0 20 atoms/cm 2 ), measured ~ 0.4 deg away 
from IGR J19405-3016. We therefore fix the column density 
at 8.73xl0 20 atoms/cm 2 and fit again the XRT data with 
a power law model. We have the fittings with the reduced 
X 2 derived around 1.1 (see Table 2). We find that the XRT 
fluxes in these snapshots can vary up to 39% with respect 
to the average at soft X-rays. The overall spectrum is rela- 
tively soft irrespective of the change in flux. The summed 
XRT data can be well-fitted by using a simple power law 
model, with a reduced x 2 ~ 1-13 for 198 do0. This suggests 
that the spectral evolution, if any, should be not very strong 
over the three XRT snapshots. See Table 2 for details of the 
spectral results. 



1 http : //heasarc . gsf c . nasa. gov/ cgi-bin/Tools/W3nh/w3nh 

2 degrees of freefom 



3.3. Combined XRT/ISGRI spectrum 

The summed XRT spectrum is combined with the ISGRI 
spectrum extracted from the all INTEGRAL observations. 
The joint spectrum can be fitted with a model of simple 
power law plus fixed absorption (reduced x 2 ^1.14 for 202 
dof) (Fig. 3). The photon index T and normalization were 
derived as 2.11±0.03 and 4.70±0.10x 10" 3 ph cm~ 2 s^ 1 
keV -1 , respectively. A constant is introduced to account 
for the difference in normalization between Swift /XRT and 
INTEGRAL /ISGRI and the sporadic snapshots of XRT ob- 
servations, and is derived as 2.19^0 60' which is slightly 
more than unity. Such deviation can be ascribed to flux 
variability at soft X-rays observed in the sporadic snap- 
shots of XRT. 

4. Discussion and summary 

As listed in the 3rd IBIS catalog (Bird et al. 2007), IGR 
J19405-3016 belongs to a group of sources with significance 
less than 5 a. About 10-20 percent of them might result 
from false detections (Bird et al. 2007). By taking the cur- 
rently available INTEGRAL observations, the exposure of 
which is roughly a factor of four more than used in the 3rd 
IBIS catalog (Bird et al. 2007), we find IGR J19405-3016 
at a significance level of ~ 9.4 a at 20-60 keV in the sum 
of the observations between March 2003 and March 2008. 
The source is consistently detected at significance level iS 4 
(j, in the three observational groups covering a time period 
of almost 5 years and in the adjacent four energy bands 
between 20-60 keV. We therefore confirm the detection of a 
real source in a previous report (Bird et al. 2007) . Hence the 
caveat mentioned in Masetti et al. (2008) in optical identifi- 
cation of IGR J19405-3016 has disappeared. A similar result 
is shown as well in Beckmann et al. (2009), where a signifi- 
cance of 11.8 cr was reported for IGR J19405-3016 in 18-60 
keV. Given the distance of a redshift z=0.052 (Masetti et 
al. 2008), the source luminosity is estimated as 2.5x10 
ergs/s (1-100 keV, use H Q =70 km/s/Mpc, q =0, fi A =0.73). 
Our analysis shows that the source was rather stable over 
5 years at energies above 20 keV. The source location and 
error circle are updated accordingly, which are more indica- 
tive of a source observed by Swift/XRT at soft X-rays as 
its real counterpart. 

At soft X-rays we find that IGR J19405-3016 has a rel- 
atively soft spectrum and the XRT flux can vary up to 39% 
with respect to the average at soft X-rays. Given the low 
sensitivity of ISGRI (about 4-5 mCrab for individual SCW) 
and low average flux level at > 20 keV (~1.5 mCrab for IGR 
J19405-3016), such a variability at hard X-rays cannot be 
studied because of the lack of sufficient signal-to-noise ra- 
tio in the ISGRI data, although short variability cannot be 
excluded. In Table 1 one sees that the source is stable on a 
long time scale at hard X-rays. It has been reported several 
times that, for Seyfert 1, their fluxes can be less variable 
at hard X-rays than at soft X-rays (Beckmann et al. 2007; 
Molina et al. 2009; Gliozzi et al. 2003). A possible explana- 
tion for this is that the large number of the scatterings of 
the soft X-rays in the corona region can let the variability 
be washed out at hard X-rays. Such an idea was proposed 
in Gliozzi et al. (2003), but so far the detailed modelings 
are not yet available. 

pi The combined XRT/ISGRI energy spectrum is well- 
fitted using a model of simple power law plus fixed absorp- 
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Table 1. IBIS/ISGRI observations log of IGR J19405- Table 2. Swift /XRT observations log of IGR J19405-3016. 
3016. 



Revs. 


Date 


Expo. 


sew 


Flux 


Sig. 




MJD 


ks 




mCrab 


a 


0056-0244 


52729-53291 


262 


146 


1.5±0.3 


4.9 


0258-0371 


53332-53672 


218 


129 


1.4±0.4 


3.5 


0416-0498 


53805-54050 


879 


275 


1.5±0.2 


7.8 


0543-0664 


54184-54547 


17 


16 


<3.3 


0.0 


0056-0664 


52729-54547 


1400 


566 


1.40±0.15 


9.4 



Note: The flux and the significance are presented in the energy 
band 20-60 keV for the data in which the source was within 
an offset angle of 10 degrees. The error bars are 1 a and the 
upper limit is 2 a. 



tion component. The spectral index falls into the range of 
the spectral index of the AGN measured by Swift /BAT at 
14-195 keV band (Tueller et al. 2008). The analysis of the 
first 9 months of the data of Swift /BAT survey of AGN 
resulted in 103 AGN being detected in 14-195 keV. The av- 
erage spectral index of these samples is 1.98, with an rms of 
0.27. For IGR J19405-3016, the spectral index of 2.11±0.03 
is consistent with the media value from BAT. About 74 
BAT AGNs have the archival spectrum at soft X-rays. A 
comparison of these samples between hard and soft X-rays 
shows that the BAT spectral slope is in general 0.23 steeper 
than in the soft X-rays (Tueller et al. 2008). A possible ex- 
planation for this could be that the spectrum is closer to 
the intrinsic one at hard X-rays than at soft X-rays due 
to the influence of the material local to AGN (Nandra et 
al. 1999; Tueller et al. 2008). The broad-band (1-100 keV) 
analysis of 36 Seyfert 1 AGNs, detected by INTEGRAL in 
the 20-40 keV band, presents a power law shape with an 
average photon index of 1.7 and a dispersion of 0.2 (Molina 
et al. 2009). The photon index of 2.11±0.03 derived for 
IGR J19405-3016 in - 1-100 keV is thus slightly more than 
the result of Molina et al. (2009). A similar discrepancy is 
found with the results of Beckmann et al. (2009) regard- 
ing Seyfert 1 AGNs. The column density as measured in 
IGR J19405-3016 is rather low, and from the X-ray spec- 
tral data, there is no need to introduce a further absorption 
local to the AGN. A possible explanation for the derived 
spectral properties of IGR J19405-3016 may be that the 
low level of the column density allows for the 'true' spec- 
trum to also be detected at soft X-rays, and thus the overall 
measured spectral shape of the source is actually the true 
one. 

In summary, the most recent analysis carried out with 
the latest software releases confirms that the previous weak 
detection, hence the low-significance source IGR J19405- 
3016 as listed in Bird et al. (2007), is not spurious. As a re- 
sult, the optical identification of a Seyfert 1.2 by Masetti et 
al. (2008) is strengthened. A spectral index ~ 2.11 supple- 
mented with a low absorption may indicate that the 'true' 
spectrum of IGR J19405-3016 extends back to soft X-rays. 
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Model Observation (date, exposure, ID) 

Par 2007-07-31 2007-08-07 2008-06-07 2007-2008 

7.6 ks 5.9 ks 1.6 ks 15.1 ks 

0036657004 0036657005 0036657006 all 

r 2.07±0.05 2.11±0.04 2.19+^3 2.11±0.03 

N 3.98±0.13 5.95±0.17 3.16±0.26 4.70±0.10 

flux l-69±ooe 2.44±0.07 1.19±0.10 1.94±0.04 

X 2 /dof 1.11/95 1.12/104 1.16/15 1.13/198 



Note: The parameters F (photon index) and N (Normalization 
at 1 keV, in units of 10 -3 ph cm -2 s _1 keV -1 ) are the 
fit results at 0.3-7 keV, using a power law model. The col- 
umn density is fixed at 8.73xl0 20 atoms/cm 2 . The x 2 /dof 
is shown as well for each fit. The flux is calculated in the 
0.3-10 keV band, in units of 10 -11 erg cm' 2 s -1 . 
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Fig. 1. The 20-60 keV significance map of IGR J19405- 
3016 as obtained from the combined revolutions 0056 - 0664 
(MJD 52729-54547). The contours start at a significance 
level of 4 a with steps of 1 a. 
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Fig. 2. Sm/t/XRTskymap of IGR J19405-3016 detected in 
the three combined XRT observations. The hot spot shows 
the soft X-ray source detected by Swift /XRT and the color 
bar shows the scale associated with the image counts. The 
open square (green) and the cross symbols (red) are the 
locations of IGR J19405-3016 estimated in Bird et al. (2007) 
and in the current work, respectively. The two circles are 
the corresponding 90% confidence level error locations. 
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data and folded model 




Energy (keV) 



Fig. 3. Upper panel: the combined XRT/ISGRI en- 
ergy spectrum of IGR J 19405-30 16 obtained from the 
INTEGRAL observations of the revolutions 0056 - 0664 
and from the three combined Swift/XRT observations (ID: 
00036657004, 00036657005, 00036657006). The data are fit- 
ted by using a power law shape plus an absorption compo- 
nent fixed at 8.73xl0 20 atoms/cm 2 . The residuals (lower 
panel) show the quality of the fit. 
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